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Abstract — The U.S. Air Force currently is in the process of 
developing an Airborne Network (AN) to provide support to 
its combat aircrafts on a mission. The reliability needed for 
continuous operation of an AN is difficult to achieve through 
completely infrastructure-less mobile ad hoc networks. In this 
paper we first propose an architecture for an AN where airborne 
networking platforms (ANPs - aircrafts, UAVs and satellites) form 
the backbone of the AN. In this architecture, the ANPs can be 
viewed as mobile base stations and the combat aircrafts on a 
mission as mobile clients. Availability of sufficient control over the 
movement pattern of the ANPs, enables the designer to develop a 
topologically stable backbone network. The combat aircrafts on 
a mission move through a space called air corridor. The goal of 
the AN design is to form a backbone network with the ANPs with 
two properties: (i) the backbone network remains connected at 
all times, even though the topology of the network changes with 
the movement of the ANPs, and (ii) the entire three dimensional 
space of the air corridor is under radio coverage at all times by 
the continuously moving ANPs. 

In addition to proposing an architecture for an AN, the con- 
tributions of the paper include, (i) development of an algorithm 
that finds the velocity and transmission range of the ANPs so that 
the dynamically changing backbone network remains connected 
at all times, (ii) development of a routing algorithm that ensures 
a connection between the source-destination node pair with the 
fewest number of path switching, (iii) given the dimensions of 
the air corridor and the radius of the coverage sphere associated 
with an ANP, development of an algorithm that finds the fewest 
number of ANPs required to provide complete coverage of the 
air corridor at all times, (iv) development of an algorithm that 
provides connected-coverage to the air corridor at all times, and 
(v) results of experimental evaluations of our algorithms, (vi) 
development of a visualization tool that depicts the movement 
patterns of the ANPs and the resulting dynamic graph and the 
coverage volume of the backbone network. 



I. Introduction 

Efforts are currently underway in the U.S. Air Force to 
utilize a heterogeneous set of physical links (RF, Optical/Laser 
and SATCOM) to interconnect a set of terrestrial, space and 
highly mobile airborne platforms (satellites, aircrafts and Un- 
manned Aerial Vehicles (ANPs)) to form an Airborne Network 
(AN). The design, development, deployment and management 
of a network where the nodes are mobile are considerably 
more complex and challenging than a network of static nodes. 

"•"This research is supported in part by a grant from the U.S. Air Force 
Office of Scientific Research under grant number FA9550-09- 1-0120. 



This is evident by the elusive promise of the Mobile Ad-Hoc 
Network (MANET) technology where despite intense research 
activity over the last fifteen years, mature solutions are yet to 
emerge [1], [2]. One major challenge in the MANET envi- 
ronment is the unpredictable movement pattern of the mobile 
nodes and its impact on the network structure. In case of an 
Airborne Network (AN), there exists considerable control over 
the movement pattern of the mobile platforms. A senior Air 
Force official can specify the controlling parameters, such as 
the location, flight path and speed of the ANPs to realize 
an AN with desired functionalities. Such control provides the 
designer with an opportunity of develop a topologically stable 
network, even when the nodes of the network are highly 
mobile. We view the AN as an infrastructure (a wireless 
mesh network) in the sky formed by mobile platforms such 
as aircrafts, satellites and UAVs to provide communication 
support to its clients such as combat aircrafts on a mission. 
Just as an Airborne Warning and Control System (AWACS) 
aircraft plays a role in a mission by providing communication 
support to fighter aircrafts directly engaged in combat, we 
believe that the aircrafts and ANPs forming the AN will 
provide similar support to the combat aircrafts over a much 
larger area. As shown in Fig. 2(a), the combat air crafts on 
a mission fly through a zone referred to as an air corridor. 
In addition to forming a connected backbone network, the 
ANPs are also required to provide complete radio coverage 
in the air corridor so that the combat aircrafts, irrespective of 
their locations within the air corridor, have access to at least 
one backbone node (i.e., an ANP) and through it, the entire 
network. Accordingly, the AN is required to have two distinct 
properties: (1) the backbone network formed by the ANPs 
must remain connected at all times, even though the topology 
of the network changes with the movement of the ANPs, and 
(2) the entire three dimensional space of the air corridor is 
covered at all times by the continuously moving ANPs. To 
the best of our knowledge this is the first paper that proposes 
an architecture for an AN and provide solutions for the all time 
connected-coverage problem of a three-dimensional space with 
mobile nodes. 



One of the pioneering results in three dimensional coverage 
problem for sensor networks was presented by Haas et al 
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[3], [4] in which they concluded that the truncated octahedron 
has the highest volumetric quotient (the ratio of the volume 
of a polyhedron to the volume of its circumsphere) among 
all the space filling polyhedrons and utilized this to develop 
placement strategies for three dimensional underwater sensor 
networks. Their scheme is a centralized one. Distributed pro- 
tocol of achieving three dimensional space coverage is found 
in the research of Tezcan et ah [5]. Poduri et ah [6] later on 
introduced the notion of NET graphs and used it to obtain 
three dimensional sensor coverage. Similar research aiming at 
the coverage problem in 3D was also presented by [7]-[9]. 
However none of these researchers put any emphasis on the 
problem of obtaining coverage while the constituting nodes 
are mobile in a three dimensional space. The mobile nature of 
the ANPs in airborne networks add yet another dimension of 
difficulty to the 3D coverage problem. 

In this paper we first propose an architecture for an AN 
where airborne networking platforms (ANPs - aircrafts, UAVs 
and satellites) form the backbone or mobile base stations of 
the AN, and the combat aircrafts on a mission function as 
mobile clients. We then proceed to determine the the number 
and initial location of the ANPs, their velocity and transmis- 
sion range, so that the dynamically changing network retains 
properties (1) and (2) mentioned in the previous paragraph. 
The rest of the paper is organized as follows. In Section II, we 
provide the system model and an architecture of an Airborne 
Network. Section III formally states the connectivity problem 
for an AN. In Section IV, we provide an algorithm that finds 
the velocity and the transmission range of the ANPs so that 
the dynamically changing network remains connected at all 
times. Section V presents a routing algorithm that ensures a 
connection between the source-destination node pair with the 
fewest number of path switching. Given the dimensions of the 
air corridor and the radius of the coverage sphere associated 
with an ANP, Section VI formulates the coverage problem for 
the air corridor. Section VII presents an algorithm that finds the 
fewest number of ANPs required to provide complete coverage 



of the air corridor at all times. The Section VIII combines 
results of Sections IV and VII and presents an algorithm to 
provide connected-coverage to the air corridor at all times. 
In Section IX we briefly describe a visualization tool that 
we developed to demonstrate the movement patterns of the 
ANPs and its impact on the resulting dynamic graph and 
the coverage volume of the backbone network. The results 
of experimental evaluations of our algorithms and related 
discussion is presented in Section X. Section XI concludes 
the paper. 

II. System Model and Architecture 

A schematic diagram of our view of an AN is shown in 
Fig. 1. In the diagram, the black aircrafts are the Airborne 
Network Platforms (ANP), the aircrafts that form the infras- 
tructure of the AN (although in Fig. 1, only aircrafts are shown 
as ANPs, the UAVs and satellites can also be considered as 
ANPs). We assume that the ANPs follow a circular flight path. 
The circular flight paths of the ANPs and their coverage area 
(shaded spheres with ANPs at the center) are also shown in 
Fig. 1. Thick dashed lines indicate the communication links 
between the ANPs. The figure also shows three fighter aircrafts 
on a mission passing through space known as air corridor, 
where network coverage is provided by ANPs 1 through 5. 
When the fighter aircrafts are at point PI on their flight path, 
they are connected to ANP4 because point PI is covered by 
ANP4 only. As the fighter aircrafts move along their flight 
trajectories, they pass through the coverage area of multiple 
ANPs and there is a smooth hand-off from one ANP to another 
when the fighter aircrafts leave the coverage area of one ANP 
and enter the coverage area of another. The fighter aircrafts are 
connected to an ANP as long as they are within the coverage 
area of that ANP At points PI, P2, P3, P4, P5 and P6 on their 
flight path in Fig. 1, the fighter aircrafts are connected to the 
ANPs (4), (2, 4), (2, 3, 4), (3), (1, 3) and (1), respectively. 

One major difference between the wireless mesh networks 
deployed in many U.S. cities [10] and the ANs is the fact that. 
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while the nodes of the wireless mesh networks deployed in the 
U.S. cities are static, the nodes of an AN are highly mobile. 
However, as noted earlier, the AN designer has considerable 
control over the movements of the mobile platforms forming 
the AN. She can decide on the locality where the aircraft/ANPs 
should fly, its altitude, flight path and speed of movement. 
Control over these four important parameters, together with 
the knowledge of the transmission range of the transceivers on 
the flying platforms, provides the designer with an opportunity 
for creating a fairly stable network, even with highly mobile 
nodes. In this paper, we make a simplifying assumptions that 
two ANPs can communicate with each other whenever the 
distance between them does not exceed the specified threshold 
(transmission range of the onboard transmitter). We are well 
aware of the fact that successful communication between 
two airborne platforms depends not only on the distance 
between them, but also on various other factors such as (i) 
the line of sight between the platforms [11], (ii) changes in 
the atmospheric channel conditions due to turbulence, clouds 
and scattering, (iii) the banking angle, the wing obstruction and 
the dead zone produced by the wake vortex of the aircraft [12] 
and (iv) Doppler effect [13] . Moreover, the transmission range 
of a link is not a constant and is impacted by various factors, 
such as transmission power, receiver sensitivity, scattering loss 
over altitude and range, path loss over propagation range, loss 
due to turbulence and the transmission aperture size [12]. 
However, the distance between the ANPS remains a very 
important parameter in determining whether communication 
between the ANPs can take place, and as the goal of this 
research is to understand the basic and fundamental issues of 
designing an AN with twin invariant properties of coverage 
and connectivity, we feel such simplifying assumptions are 
necessary and justified. Once the fundamental issues of the 
problem are well understood, factors (i) through (iv) can be 
incorporated into the model to obtain a more accurate solution. 

HI. Design for Connectivity - Problem 
Formulation 

It is conceivable that even if the network topology changes 
due to movement of the nodes, some underlying structural 
properties of the network may still remain invariant. A 
structural property of prime interest in this context is the 
connectivity of the dynamic graph formed by the ANPs. We 
want the ANPs to fly in such a way, that even though the 
links between them are established and disestablished over 
time, the underlying graph remains connected at all times. 
Although we give connectedness of the graph as an example 
of a structural property, many other graph theoretic properties 
P can be specified as design requirements for the network. 
The problem can be described formally in the following way. 

Consider n nodes (flying platforms) in an m-dimensional 
space (for ANP network scenario m = 3). We denote by 
Xi{t) G the coordinates of the node i at time t, where 
by convention Xi is considered a m x 1 column vector, and 
by x(t) = [xi^(t), . . . , x^^(t)] , the mn vector resulting 
from stacking the coordinates of the nodes in a single vector. 
Suppose that the dynamics of node i (for alH G {1,2,..., n}). 



is given by Xi{t) = Ui{t), where Ui{t) is the control vector 
taking values in some set U C W^. In vector notation, the 
system dynamics become 

x(t) = u(t) (1) 

where x(t) = [xf (t), . . . , xj(t)]^ and u(t) = 
[ui^{t)^ . . . ^Un^{t)]^ are mn x 1 vectors, respectively. The 
network of flying platforms described by system (1), gives rise 
to a dynamic graph Q{^{t)). 

Q{^{t)) = (V,f (x(t))) is a dynamic graph consisting of 

• a set of nodes V = {1, 2, . . . , n} indexed by the set of 
flying platforms, and 

• a set of edges f (x(t)) = {(i, j) | dij{^{t)) < 6} with 
dij{-K{t)) =11 Xi{t) — Xj{t) II as the Euclidean distance 
between the platforms i and j and ^ > is a constant. 

Since we have control over the node dynamics, the question 
that naturally arises is whether we can control the motion of 
the ANPs so that ^(x(t)) retains graph-theoretic properties 
of interest P for all time ^ > 0. A graph Q is connected 
if there exists a path between any two nodes of the graph Q. 
Often times the property P will correspond to the requirement 
that the graph Q remains connected at all times. Formally the 
problem can be stated as follows. Suppose that Cn,p is the 
set of all graphs on n nodes with property P. Is it possible 
to find a control law u{t) such that if 5(x(0)) G Cn,p then 
a(x(t)) G Cn,P for aU t > 01 

Although a few researchers have studied problems in this 
domain [14]-[16], many important questions still remain unan- 
swered. For example, in our study of the movement pattern 
of the ANPs to create a connected network, we assume that 
the flight paths of the mobile platforms are already known and 
we want to find out the speed at which these platforms should 
move, so that the resulting dynamic graph remains connected 
at all times. The studies undertaken in [14]-[16] do not address 
such issues. Although the movement of the airborne platforms 
will be in a three dimensional space, in a simplified version 
of the problem in two dimension (i.e., when all the aircrafts 
are flying at the same altitude) the problem can be stated as 
follows: 

Mobility Pattern for Connected Dynamic Graph (MPCDG): 
This problem has five controlling parameters: 

(i) a set of points {pi,^>2, • • • ,Pn} on a two (or three) 
dimensional space (representing the centers of circular flight 
paths of the platforms), 

(ii) a set of radii {ri, r2, . . . , r^} representing the radii of 
circular flight paths, 

(iii) a set of points {/i, ^2, • • • , ^n} representing the initial 
locations (i.e., locations at time t = 0) of the platforms on 
the circular flight paths, 

(iv) a set of velocities {t'l, . . . , t'n} representing the speeds 
of the platforms, and 

(v) transmission range of the transceivers on the airborne 
platforms. 

IV. Design for Connectivity - Solution 

In the MPCDG problem scenario, any structural property 
P of the resulting dynamic graph will be determined by the 
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Fig. 3. Vector representations (Ri(t) and Rj(t)) of two points i and j at 
time t moving along two circular orbits: = 15, rcj = 27, ZciOx = 



problem parameters (i) through (v). The problems that arise 
in this formulation are as follows: Given any four of the five 
problem parameters, how to determine the fifth one, so that 
the resulting dynamic graph retains property P at all times? 
Most often we would like to know that given (i), (ii), (iii) and 
(v), at what speed the ANPs should fly so that the resulting 
graph is connected at all times. Alternately, we may want to 
determine the minimum transmission range of the ANPs to 
ensure connectivity. In this case, the problem will be specified 
in the following way. Given (i), (ii), (iii) and (iv), what is the 
minimum transmission range of the ANPs so that the resulting 
graph is connected at all times? In order to answer these 
questions, we first need to able to answer a simpler question. 
Given all five problem parameters including the speed of the 
ANPs, how do you determine if the resulting dynamic graph 
is connected at all times? We discuss this problem next. 

A. Connectivity checking for AN when all control parameters 
specified 

In this subsection we describe our technique to find answer 
to the question posed in the previous paragraph. Suppose 
that two ANPs, represented by two points i and j (either in 
two or in three dimensional space, the two dimensional case 
corresponds to the scenario where the ANPs are flying at same 
altitute) are moving along two circular orbits with centers at 
Ci and Cj with orbit radius Vi and Vj as shown in Fig. 3 with 
velocities Vi and Vj (with corresponding angular velocities uji 
and ujj), respectively. 

A moving node i is specified by the radius vector Ri{t) 
directed from some origin point O, and similarly Rj[t) for 
point j. Therefore the distance Sij{t) between the nodes i— j 
at time t is given by: 



slit) 



(Riit) - Rj{t)f = Rl{t) + R^t) - 2Ri{t) ■ Rj{t) 

(2) 

As mentioned earlier, we have assumed that the communica- 
tion between the ANPS is possible if and only if the Euclidean 
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Fig. 4. Initial phase angle of point i\ at time point is shown as 2(0) 

distance between them does not exceed the communication 
threshold distance D. This implies that the link between the 
nodes i and j is alive (or active) when 



(3) 



In the analysis that follows, we have assumed that ANPs 
are flying at the same altitude, i.e., we focus our attention 
to the two dimensional scenario. However, this analysis can 
easily be extended to the three dimensional case to model the 
scenario where the ANPs are flying at different altitude. In this 
case we can view the ANPs as points on a two-dimensional 
plane moving along two circular orbits, as shown in Fig. 3. 
In Fig. 3, the vectors from the origin O to the centers of 
the orbits q and cj are given as and Tc^.. The cartesian 
co-ordinates of the centers can be readily obtained as = 
{rc,cos ac,,rc,sin a^) and = {vcjCos ac^.Tc^sin a^). 
Accordingly, Ri{t) can be expressed in polar coordinates: 
Ri{t)^Oi{t) with respect to origin point O, as shown in Fig. 3, 
and similarly for Rj (t). The initial location of the points Ri{0) 
and Rj{0) are given. From Fig. 4, the phase angle pi for node 
i with respect to the center of orbit c^, can be calculated as 
(by taking projection on the axes): 

Ri{0)cos Oi{0) 



tan j3i = 
Since from Fig. 3, 



Ri{0)sin Oi{0) — Tc^sin 



Ri{t) 



(4) 



(5) 



where fi(t) = (r^ cos (/^^ + cj^t), sin (/^^ + cj^t)) (since 
angle made by i at time t w.r.t. q is given by {Pi -\- ujit)). 
Therefore, the angle between r^(t) and Tc- is {Pi — a^^ ^-ooit). 
Hence, 

R^if) = rl + rf + 2rcji cos (A - a^, + cuit) (6) 

Now taking the projection of Ri{t) = fa -^u{t) on the x and 
y axes, we get 

Ri{t) cos Oi{t) = cos + cos (ft + cj^t), (7) 
Ri{t)s\ii Oi{t) = Tc- sin ac- + sin (ft + cj^t) (8) 

Recalling cos (A — B) = cos ^4 cos 5 + sin A sin 5, and 
simplifying we get 

Ri{t)Rj{t) cos{Oi{t) - 6j{t)) = Vcjcj cos a^cj 

-\-rirj cos(ftj + {coi - ujj)t) + r^r^ cos(cec^ - Pj - ^jt) 

-{-rcjVi cos{acj - Pi- ujit) (9) 
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where Ofc,,- = ac^-ac^ and Pij = Pi- 
2 with equations 6 and 9, we have: 



-Pj. Combining equation 



+ ^c,- + + 2rc^.r_^- cos(/3_^- - a^- + ^j^) 

+ rc,rc^. cos ac^cj + cos(/3ij + {uji - Uj)t) 

+ rciVj cos (c^ci ~ y^j ~ ^jt) 

+ Tc^.r^ cos(ac^. - A - ^it) (10) 

In equation 10, all parameters on the right hand side are 
known from the initial state of the system, and thus the 
distance Sij (t) between the nodes i — j at any time t can be 
obtained. If the ANPs move at same velocity, i.e., cj^ = ooj = cu 
for all j and the radius of the circular orbits are identical, i.e., 
ri = Tj = r for all i, j, and the above expression simplifies 



to: 



^iji'^) = + + 2rc.r cos(/3i — ac- + ujt) 

+ Vcjcj cos a^c- + cos ft J 

+ rc,rcos{ac, - - cot) 

+ r cos(ac,- - Pi- <^t) 



(11) 



If the problem parameters (i) through (v) are specified, we 
can check if the dynamic graph is connected at all times 
following these two steps. In the first step, we determine the 
lifetime (active/inactive) of a link between a pair of nodes i 
and j in the following way. 

Algorithm 1: Link Lifetime Computation 

1. begin 

2. Using equation (10), compute and plot the distance 

between a pair nodes i and j, Sij{t), as a function 
of time; ( See Fig. 5(a)) 

3. Draw a horizontal line in the Sij{t) versus t plot with 

Sij (t) = D, where D is the communication threshold, 
i.e., communication between i and j is possible 
if Sij{t) < D and impossible otherwise. 
Call this line communication threshold line, CTL. 

4. The CTL is divided in to segments corresponding 

to the parts where Sij{t) < D and where Sij{t) > D. 

5. Projections of the CTL segments on the x-axis 

(i.e., the time line) indicates the times when 
the link between i and j is alive and 
when it is not. ( See Fig. 5(b)) 

6. end 

Using Algorithm 1, we can compute the life time of every 
link (i.e., every pair of nodes) in the network. In step 2, using 
Algorithm 2 (given below) we divide the time line into smaller 
intervals and determine exactly the links that are active during 
each of these interval. For each of the intervals we check if the 
AN graph is connected during that interval using connectivity 
checking algorithm in [17]. The algorithm is described in 
detail next. 



Algorithm 2: Checking Connectivity of Airborne Network 
between time t = ti and t = t2 

1. begin 

2. Using the Algorithm for Link Lifetime Computation, 

compute the lifetimes of links between all node pairs 
and plot them over time line. (See Fig. 6) 

3. Draw a vertical line through start and finish 

time of each interval associated with a link on the 
X-axis (time line) 

4. Repeat step 3 for each link of the network 

5. The X-axis (time line) is divided into a number of 

smaller intervals. (See Fig. 6, intervals are 
numbered from 1 through 17). From the figure, 
we can identify all the links that are alive 
during any one interval. 

6. Check if the AN graph is connected with the set 

of live links during one interval. This can be done 
with the connectivity testing algorithm in [17] 

7. Repeat step 6 for all the intervals between t = ti 

and t = ^2- 

8. If the AN graph remains connected for all intervals, 

conclude that the AN remains connected during 
the entire duration between t = ti and t = t2, 
otherwise conclude that the specified problem 
parameters does not ensure a AN that remains 
connected during the entire time interval 
between t = ti and t = t2. 

6. end 

An example of a plot of equation (10) (generated using 
MATLAB) is shown in Fig. 5(a) with communication thresh- 
old distance = 18. This implies that the link between the 
nodes i and j exists, when the distance between them is at 
most 18 and the link does not exist otherwise. This is shown 
in Fig. 5(b). The red part indicates the time interval when the 
link is inactiveiov dead) and the blue part indicates when it is 
active (or live). 

Thus using equation (10) and comparing the distance be- 
tween any two nodes with the threshold distance D, we 
can determine active/inactive times of all links. This can be 
represented as intervals on a time line as shown in Fig. 6. By 
drawing projections from the end-points of the active/inactive 
times of each link on the time line, we can find out all the 
links that are active during a interval on time line. As shown 
in Fig. 6, links 1, 2 and 3 are active in interval 1; links 1 and 
3 are active in interval 2, links 1, 2 and 3 are active in interval 
3 and so on. Once we know all the links that are active during 
a time interval, we can determine if the graph is connected 
during that interval using any algorithm for computing graph 
connectivity [18]. By checking if the graph is connected at all 
intervals, we can determine if the graph is connected at all 
times, when the ANPs are moving at specified velocities. 

B. Finding the velocity of the ANPs to ensure a connected AN 
during operational time between t = ti and t = t2 

In subsection A, we have described a technique to determine 
if the AN remains connected during the entire operational time 
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(a) Distance between two points i and j as a function of time 
Fig. 5. Effect of the distance between nodes on the existence of the communication Hnk between them 



(b) Active (Blue)/Inactive (Red) times of the Hnk between i and j with 
transmission range Tr = 18 



between t = ti and t = ^2), when all problem parameters, 
(i)-(v) are specified. In this subsection we try to determine 
the problem parameter (iv) (i.e., velocity of the ANPs) that 
we ensure a connected AN during the entire operational time 
between t = ti and t = ^2) when all other problem parameters 
have already been specified. The minimum and maximum 
operating velocities of the ANPs (Vmin, ^max) are known. By 
conducting a binary search on this range, we can compute 
the minimum velocity at which the ANPs should fly, so that 
the AN remains connected during the entire operational time. 
Alternately, we can also try to determine the velocity at which 
the ANPs should fly, so that the AN remains connected during 
the entire operational time and fuel consumption by the ANPs 
is minimized. If it is known that the fuel consumption is 
minimized when the ANPs fly with velocity VopUmah we can 
find the velocity that is closest to VopUmai and also ensures 
connectivity of the AN during entire operational time by a 
targeted search within the range (Vmin^'^max)- 

C. Finding the transmission range of the ANPs to ensure a 
connected AN during operational time between t = ti and 

t = t2 

In this subsection we try to determine the problem parameter 
(v) (i.e., transmission range of the ANPs) that we ensure a 
connected AN during the entire operational time between t = 
ti and t = t2) when all other problem parameters have already 
been specified. The maximum transmission range of an ANP 
is known in advance (Tmax- By conducting a binary search 
within the range — Tmax^ we can determine the the smallest 
transmission range that will ensure a connected AN during 
the entire operational time when all other problem parameters 
have already been determined. 

V. Routing with Minimum Path Switching 

In the previous section we described a procedure to deter- 
mine the velocity of the ANPs so that the resulting dynamic 
graph is connected at all times. Although the graph remains 
connected at all times, as the links come and go (alive or dead) 
a path between a source-destination pair may not exist for the 
entire duration of communication. Suppose that a node s has to 
communicate with another node d from time t = 5tot = 12. 
Since the graph is connected at all times, at least a path, say 
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Fig. 6. Active/Inactive time interval of each link and interval intersection 
projections on the time line 



Pi, exists from s io d dX t = 5. However, this path may 
not exist till t = 12. Suppose that as one of its link dies. Pi 
breaks at t = 7. Clearly Pi cannot be used for communication 
between s and d dXt = 7. Since the graph is connected at all 
times, there must exist at least one path, say P2, between s 
and dsitt = 7. Therefore data can be transferred from s to d 
using P2 at t = 7. However, P2 can break at t = 10, in which 
case a yet another path, say P3 (which is guaranteed to exist 
because the graph is connected at all times) can be used for 
communication between s and d from t = 10tot = 12. In 
such a scenario, the path sequence Pi ^ P2 ^ P3 is used 
for communication between s and d in the time interval t = 5 
to t = 12. In this scenario the path has to be switched two 
times, once from Pi P2 and the other time from P2 ^ P3. 
However, it is possible that communication between s and d in 
the time interval t = 5 to t = 12 could have been achieved by 
only one path switching, using a path P4 from t = b to t = 9 
and a path P5 from t = 9 to t = 12. Since path switching 
involves a certain amount of overhead, it is undesirable and 
as such we would like to accomplish routing for the duration 
of communication with as few path switching as possible. 

In Fig. 6 we showed how "lifetime" of a link (i.e., 
alive/dead) can be computed. Since paths comprise of links, 
a path between a source-destination node pair will also be 
alive/dead at different points of time. From the lifetime of 
links, we can compute the lifetime of paths in the following 
way. If the number of nodes (i.e., ANPs) in the network is n, 
there exists n(n — l)/2 links with each having an individual 
lifetime. If a path P is made up of links /i, ^2, • • • , ^fe, the 
path P is "alive" when all the links h through 1^ are alive. 
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Fig. 7. Lifetime of r Paths between a source-destination pair, and the pig. 8. Solution Produced by the Optimal Algorithm and the Minimum Path 
corresponding time intervals when they are alive, s and d need to communicate Switch Routing Algorithm 
between time t = ti and t = t2 



Therefore, similar to Fig. 6 that shows the"Hfetime" of a Hnk, 
we can construct a figure for 'lifetime" of a path (Fig. 7). Once 
we have knowledge of lifetimes of paths, we can construct a 
route from the source node s to the destination node d with 
the fewest number of path switching in the following way. 

The lifetimes of paths between a source-destination node 
pair s and d is shown in Fig. 7. The time intervals 
during which a path is alive is shown by solid lines in 
Fig. 7. We use the notation Pj = {/j^i, /j^2, • • • , 
to indicate that the path Pj is alive during the jk time 
intervals {/j^i, /j^2, • • • , ^jj^}' shown in the Fig. 7. In 
the application scenario that we are considering, we want a 
communication channel to be open between the source node 
s and the destination node d for the entire duration of time 
from t = ti to t = t2. Since it is possible that no single path 
between sto d remains alive for the entire duration from t = ti 
to t = ^2, a set of paths V may constitute a communication 
channel from s to d for the duration, where each path in V 
is alive only for a fraction of the time interval from t = ti to 

t = t2. 

Next we focus on the number of paths between s and d 
that we need to consider. Since the graph has n nodes and 
n{n — l)/2 links, there could be as many as (1 + (n — 2) + 
(n - 2)(n - 3) + . . . + (n - 2)(n - 3)(n - 4) ... 2 = 0(n^)) 
paths corresponding to 1-hop, 2-hop, . . ., (n — l)-hop paths 
between s and d. It may also be noted that each of these 
paths will have a lifetime associated with it. Since examining 
O(n^) paths and their lifetimes will be too time consuming, 
we restrict our attention to only those paths between s and d 
whose number of hops is at most k, for some specified value 
of k. Restricting the number of hops in a path to at most k 
(from n — 1), we reduce the computational complexity from 
O(n^) to 0{n^). Suppose the set of paths of at most k hops 
is denoted at Vk- The minimum path switch routing algorithm 
given below finds a subset C Vj^ so that paths in Vj^ 
maintain a communication channel between s and d for the 
entire time duration from t = ti to t = t2 with the fewest 
number of path switchings. 

Minimum Path-Switch Routing Algorithm 

Input: The set Vk of paths between s and d with length at most 
/c-hops, and associated lifetimes of each each path Pi e Vk, 
in the form of live intervals of Pi. If there are r live intervals 
of Pi in the time interval between t = ti to t = ti, we denote 



Pi — h,2-) • • • •) h,r}- 

Output: A subset C Vk so that paths in V'^ maintain a 
communication channel between s and d for the entire time 
duration from t = ti to t = t2 with the fewest number of path 
switchings. 

Comments: The algorithm uses a greedy (locally optimum) 
approach to find the paths needed to have one live path from 
s to d during the entire time interval t = ti and t = t2. 
In Theorem 1, we prove that this locally optimum greedy 
approach indeed finds the globally optimal solution. 

1. begin 

2. 7^^ = 0; 

3- t start ~ 

4. ^finish ^1? 

5. While (t finish < t2) do 

begin 

(i) for all Pi e Vk do 

((start_time {lij) < tstart) 
&& {t finish < finish_time{Iij) for some j) 
begin 

(i) t finish = finish_time{Iij) 

(ii) vj, = r,uPi- 

end 

tstart ~ tfinish") 

end 

6. end 

Theorem 1: The Minimum Path-Switch Routing Algorithm 
finds a set of paths so that a communication channel is open 
during the entire duration from t = ti to t = ti with the 
fewest number of path switches. 

Proof: Suppose that an optimal algorithm selected the paths 
{Pi^^Pi^^ . . . , Pi^} and the Minimum Path-Switch Routingl 
algorithm selected the paths {P^/^ , P^^ , . . . , P^/^ }, where x < y. 
The live intervals of the paths {Pi^^Pi^^ . . . ^ P^^ } that were 
selected by the optimal algorithm are {U^ , /i2,i2 ^ • • • ^ } 
respectively. Similarly, the live intervals of the paths 
{P^/^ , Pi>^, . . . , Pi>^} that were selected by the Minimum Path- 
Switch Routing algorithm are {/^/^ , /^/^ . . . , 7^/^^^/^} re- 
spectively. The paths and intervals chosen by the two al- 
gorithms are shown in Fig. 8. As shown in Fig. 8, the 
finish times of the intervals {Iii,jx^Ii2,j2^ • • • ^^ixdx} 
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denoted as {Ti^j^^Ti^j^^ . . . ^Ti^j^} and the finish times 
of the intervals {/^/^ , /^^ j/ , . . . , /^/^ j^} are denoted as 

Since the Minimum Path-Switch Routing Algorithm 
chooses the path Pi G Vk such that Pi is live at tgtart and the 
finish time of the live interval containing is largest among 
the finish times of the intervals associated with all the paths, 
we can conclude that Ti'^j[ > Ti^ji- Therefore, replacing the 
path Pi^ from the optimal soultion by the path P^/^ we will have 
a new optimal solution {P^/^ , P^2 , . . . , Pi^}- Because of nature 
of the path selection criteria of the Minimum Path-Switch 
Routing Algorithm, we can conclude that > ^22,^2- 

Therefore, replacing the path Pi^ from the new optimal 
soultion by the path P^^ we will have yet another optimal 
solution {P^/^ , P^^ , . . . , Pi^ }. Continuing this process, we can 
get an optimal solution {Pi>^^Pi^^ . . . , Pi^}. This implies that 
the Minimum Path-Switch Routing Algorithm will select only 
X paths instead of y, (x < y) to have an open communication 
channel for the enire duration of t = ti to t = t2, and since 
X is the optimal number of paths needed for this purpose, 
the Minimum Path-Switch Routing Algorithm produces an 
optimal solution. 

VI. Design for Coverage - Problem Formulation 

In this section, we discuss the coverage model of the 
network formed by the ANPs. As shown in Fig. 2(a), an 
air corridor through which the combat aircrafts fly towards 
their destination can be modeled as a collection of rectangular 
parallelepipeds. As the combat aircrafts must have access to 
the AN as they fly through the air corridor, all points inside 
the air-corridor must have radio coverage at all times. As 
the shape of an air corridor can be quite complex, we view 
that the complex shape can be approximated with rectangular 
parallelepipeds as shown in Fig. 2(b). 

The length, width and height of this section of the air- 
corridor are denoted as Lac, Wac and Hac, respectively. The 
radius of the circular orbit of the ANPs is denoted by Tq and 
the number of ANPs in each such orbit is denoted by n. We 
assume that the ANPs move around in their orbit with uniform 
velocities. The coverage volume of each ANP is defined as 
a spherical volume of radius with the ANP being at the 
center of the sphere. We assume that the orbits of the ANPs are 
located at the top surface of the air corridor so that they do not 
cause any hindrance in the flight path of the combat aircrafts. 
This is shown in Fig. 2(c), where two circular orbits, each of 
them containing 5 ANPs, are located at the top surface of the 
air corridor. The number of orbits located at the top surface 
of the air corridor section is denoted as m, and accordingly 
the total number of ANPs in the network is given by mn. 
The goal of the coverage problem is to provide complete 
coverage at all times of the entire air corridor with the fewest 
number of ANPs. In this problem, a complete coverage must 
be provided irrespective of the locations of the ANPs as they 
move continuously in their respective orbits. 

In Figs. 9, 10, 11, we can see the 3D view, top-view and 
front- view of 5 ANPs in circular orbits and the volume covered 
by them for three different cases. In case I (Fig. 9), the value of 



orbit radius (Vo) is greater than that of the spherical coverage 
volume (rs) of each ANP. It can be clearly seen from Fig. 9(b), 
that there is an open space inside the orbit, that is not covered 
by any of the ANPs as they move in the orbit. To increase the 
intersection volume, can be at most (case II). In Fig. 10, 
To = rs, and we can see from Fig. 10(b), that the spheres meet 
at a single point inside the orbit. We call the intersection of 
adjacent spheres as leaf, the top view of which is visible in 
the figure. To increase the intersection of the spheres further, 
we need to decrease ro even more, and thus bringing in the 
ANPs even more closer to each other. This is case III (shown 
in Fig. 11), where < r^. 

The coverage problem can be stated as follows: Given the 
rectangular parallelepipeds in terms of Lac, Wac and Hac and 
the radius of the coverage sphere associated with an ANP, r^, 
find the radius of the orbit of the ANPs r^, and the number of 
ANPs in each orbit (n), the entire volume of the air corridor 
is covered at all times with the fewest number of ANPs. 

Intersection of coverage spheres of the ANPs create a 
coverage volume. For two intersecting spheres, the intersection 
volume is shown in Fig. 12. As the ANPs move in their 
orbits, the associated coverage spheres move with them and 
consequently the volume that is covered by the moving ANPs 
also changes. As a consequence some volume will be covered 
only a part of the time. However, a part of the intersection 
volume will be covered at all times irrespective of the positions 
of the ANPs as they move in their orbits. This is defined as the 
invariant coverage volume. We would like to use this invariant 
coverage volume as building blocks in order to fill up the air 
corridor modeled in the form of a rectangular parallelopiped. 
Since the invariant coverage volume is irregular- shaped, it is 
difficult to use it as a building block. For ease of coverage 
using a building block with a regular shape, we extract a 
cylindrical volume out of this invariant volume and use it to 
fill up the rectangular parallelopiped. Such a cylinder is shown 
in Fig. 12. Different views of such a cylindrical section for 5 
intersecting ANPs in a circular orbit are shown in Fig. 13. 

As we have decided to use a cylinder as the building block 
to cover the air corridor, we need to know the height and radius 
of the circular surface of such cylindrical blocks, denoted by 
2hc and rc, respectively. The height and radius of the invariant 
coverage cylinder are determined by (i) the orbit radius (ro), 
(ii) the number of ANPs per orbit (n) and (iii) the radius of the 
spherical coverage volume of each ANP (r^). As mentioned 
earlier, in this design the ANPs and their orbits are placed on 
the top surface of the air corridor. As a consequence, the top 
half of the invariant coverage cylinder cannot be utilized and 
only the bottom half of the cylindrical volume (of height he) 
will be used for the coverage of the rectangular parallelopiped. 
Therefore, in order to cover the height of air corridor, one must 
satisfy the constraint he > Hac (Fig- 14). Once this constraint 
is satisfied, the problem reduces to cover the plane defined by 
Lac X Wac with circlcs of radius Tc with a goal to minimize 
the total number of ANPs required (mn). 

We investigate the structure of the coverage volume. In 
Fig. 15, we have shown the top view through the center 
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(a) 5 ANPs moving in a circular orbit (b) Top View of the coverage spheres and cir- (c) Front View of the coverage spheres 

cular orbit of 5 ANPs and circular orbit of 5 ANPs 



Fig. 9. Case I: orbit radius of ANPs (to) > radius of coverage sphere (r^) 




(a) 5 ANPs moving in a circular orbit (b) Top View of the coverage spheres and cir- (c) Front View of the coverage spheres 

cular orbit of 5 ANPs and circular orbit of 5 ANPs 



Fig. 10. Case II: orbit radius of ANPs (ro) = radius of coverage sphere (rg) 




(a) 5 ANPs moving in a circular orbit (b) Top View of the coverage spheres and cir- (c) Front View of the coverage spheres 

cular orbit of 5 ANPs and circular orbit of 5 ANPs 



Fig. 11. Case II: orbit radius of ANPs (ro) < radius of coverage sphere (rg) 



of orbit of three consecutive spheres (out of total n of 
them) intersecting with each other and moving around in the 
circular orbit with center at O. The center of the spheres are 
denoted as Ci, C2 and C3. The radius of the circular orbit is 
OCi = OC2 = OCs = ro. All the spheres are of uniform 
radius , and moving in a uniform velocity. The term leaf is 
used to refer to the intersection between two adjacent spheres. 



In the^top view, i^ can be seen as the intersection of two circular 
arcs PSQ and QRP. The length of the leaf PQ is denoted 
by 2hi (Ft = TQ = hi). Distance from the center of the 
orbit center O to the end point of the leaf Q is denoted by y. 
Therefore, 

TO = TQ-OQ = hi-y 
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Fig. 12. Intersection volume of two spheres; cylindrical 
volume cut from the intersection volume 



Cylindrical intersection 
region of spheres 




front view of cylindrical intersection region 




(a) 3-D view 



(b) Front View 



Fig. 13. Cylindrical volume cut from the intersection of the coverage spheres of ANPs 

(rs > To) 













/ 
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Fig. 14. Air-corridor being filled up with cylindrical sections - the cylinder 
shown is formed due to the ANPs in a particular orbit on the top surface 



ANPs 




\ 

Coverage 
'Sphere of ANP 



The width of the leaf (TS) is denoted by wi . Angle between 
two adjacent leaves ZPOM is given by 0. For n number of 
spheres moving in the orbit, is given by Therefore, 

ZPOC2 = ZC2OM =i = - 
2 n 

Now, in Fig. 16, the side view of two intersecting spheres 
and leafs are shown. The intersecting spheres are shown using 
dashed line, whereas the intersecting leaves are shown using 
solid lines. We cut the largest cylindrical volume from the 
intersecting region such that this is covered by at least one 
sphere at all times, as the spheres move around in the orbit. 
From Fig. 15: 

i.e.,h^ = rl-{rs-wif 

hi = ^wi{2rs - wi) (12) 

Also, from AC2TO: 



C2T = TOtan- 

TT 

Ts-wi = {hi - y)tan- (13) 

n 

0^ = C^sin^ 


Vs-wi = ToSin- 

7T 

wi = rs — roSin— (14) 
n 

From equations 12, 13 and 14, it is clear that all variables 
hi, wi and y can be represented through variables Tq and n 
only. Now from Fig. 16, it is clear that the cylindrical volume 
with largest height that remains covered during the movement 



ANPa 



Fig. 15. Top View of Three Consecutive Intersecting Spheres Moving in a 
Orbit 




Fig. 16. Side View of Intersecting Spheres and the Leaves 



of the spheres, has radius of the circular surface (rc) equal to 
length VA. This can be calculated as : 

VA = OU 

= 2{TQ-OQ) 
Le.,rc = 2{hi-y) (15) 

The height of the cylinder is given by VW = 2VO = 2hc, 
calculated as: 

hi = hf-{hi-yf 

he = ^y{2hi - y) (16) 

Therefore, given the radius of the orbit r^, the number of 
spheres moving in each orbit n, and radius of each sphere r^, 
the cylinder can be determined by using the above equations 
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number of spheres per orbit (n) 




orbit radius (r ) 



number of splieres per orbit (n) 




number of spheres per orbit (n) 



(a) Radius of Cylinder (b) Height of Cylinder (c) Volume of Cylinder 

Fig. 17. Change of Cylinder Radius, Height and Volume with Orbit Radius and Number of Spheres per Orbit with fixed Sphere Radius =10 units 



as follows: 



VII. Design For Coverage - Solution 



2{hi - y) 
2{rs - wi)cot- 



2roSin-cot- 

zrocos- - zrocos — 
2 n 



(from equation 15) 
(from equation 13) 

(from equation 14) 



(17) 



Using equations 12, 13, 14 and 16, we can express he as : 
hi = hf-{hi-yf 



= (r^ - {ts - wif) - {{ts - wi)cot-f 



2 2 



The volume of cylinder is given by: 



(18) 



(19) 



The change of the radius (rc), height (he) and volume (Vc) 
of the cylindrical region with the values of and n for fixed 
value of can be seen in the plots in Figs. 17(a), 17(b) and 
17(c). The coverage problem for an airborne network can be 
formally defined as follows. Given: 

• The length, width and height of the air corridor as Lac, 
Wac and Hac- 

• The radius of each coverage sphere associated with each 
ANP Ts. 

Find (i) the orbit radius, r^, (ii) the number of ANPs in each 
orbit, n, (iii) the number of orbits m required to cover the 
air corridor, and (iv) the placement of the center of the orbits 

{xi^Hi)^ 1 < i < m, such that: 

1) The orbits are placed only on the top surface of the air 
corridor. 

2) All points in the rectangular area defined by Lac and Wac 
is covered by at least one circle of radius rc with center 

at {xi.yi), 1 <i <m. 

3) he > Hac 

4) Total number of spheres required = total number of orbits 
X number of spheres per orbit = mn is minimized. 



In order to minimize the objective function subject to the 
constraints, we need to find Tq and n. This will determine the 
values of Tc and he. It can be seen from equations (17) and 
(18), that decreasing increases h^ but decreases Tc. Intu- 
itively, with smaller value of we will need more orbits (m) 
to cover the rectangular parallelopiped, which will eventually 
increase the number of ANPs required (mn). Therefore, for 
given Ts we need to set at the highest possible value that still 
satisfies the constraint he > Hac- For the corresponding m 
will be determined by the placement strategy of the orbits on 
the top surface of the air corridor. The overall objective of the 
coverage problem is to minimize mn subject to the constraint 
he > Hac- We use the following two strategies for placement 
of circular orbits on the top surface of the air corridor: 
Strategy 1: The largest square that can be inscribed in the 
circular surface of the cylinder is used as the building block 
to cover the rectangular region defined by Lac and Wac- With 
Tc being the radius of the cylinder, the length of each side a of 
the square is V2re- Therefore, total number of orbits can be 
calculated as m = [^^1 x T^^J^l- Hence, the optimization 
problem following strategy 1 can fee formally stated as: 



minimize mn 



La 



\ 



subject to : h, 



V2r. 



Wac 



cos- 



X n 



] X n 



-I > Ha 



(20) 
(21) 



Strategy 2: The largest rectangle that can be inscribed in the 
circular surface of the cylinder, and that has the same length 
to width ratio as that of the rectangular region defined by Lac 
and Wac, is used as the building block for covering the entire 
region. Let a and b be the length and width of such a building 
block. Tc being the radius of the cylinder, we get the following: 



b ~ Wa. 



and (2rc)^ = + 6^ 



From the above two relations we can get: 



^VcLa 



2rrWa 



Since Tc = 2roCos^, total number of orbits (cylinders) 
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(a) Placement using Strategy 1 (b) Placement using Strategy 2 

Fig. 18. Covering Lac x Wac plane with circles of radius Tc using strategy 
1 and 2 

required m to cover the entire region can be calculated as: 

Accordingly, the optimization problem following placement 
strategy 2 can be formally stated as: 

minimize mn \— — — ^r-^l^ x n (22) 

^ n 

subject to: he = V^s -^o ^ (23) 

The diagrams for the placement of orbits, and hence the cylin- 
drical regions following the above two strategies, are shown 
in Fig. 18(a) and 18(b), respectively. Hence, the location of he 
center of the orbits (xi^i/i) can be easily determined. It may 
be observed both placement strategy 1 and 2 formulates the 
coverage problem as a non-linear optimization problem. We 
used the non-linear constrained program solver Nimbus [19] 
to solve optimization problems following strategies 1 and 2. 
The results obtained from Nimbus is discussed in Section X. 

VIII. Design for Connected Coverage 

In Section VI, we discussed the three dimensional air corri- 
dor coverage problem with the ANPs. As the ANPs are mobile, 
the coverage volume associated one ANP is continuously 
changing with time. In Section VII we described techniques 
to find least cost solution to the air corridor coverage problem 
with mobile nodes. Earlier in Section IV, we discussed how 
to determine the velocity and subsequently the transmission 
range of the of the ANPs, so that resulting backbone network 
formed by the ANPs remain connected at all times. In this 
section, we discuss how to design a network of ANPs, so that 
(i) it remains connected at all times and (ii) it provides 100% 
coverage to the air corridor at all times. 

We provide a two phase solution to the connected coverage 
problem. In the first phase, using the techniques described 
in Section VII, we determine the number and orbit of the 
ANPs that will provide 100% coverage to air corridor at all 
times. Once that is accomplished, in the second phase, using 
the techniques described in Section IV, we determine the the 
velocity and the transmission range of the ANPs so that the 
backbone network formed by the ANPS remains connected at 
all times. 




Fig. 19. Snapshot of Visualization Tool with Three Moving Objects 

IX. Visualization Tool for Airborne Network 
Design 

The visualization tool was designed for observing the 
movement of objects along circular orbits in a 3D plane. 
The Euclidean distance between every pair of objects keeps 
changing due the their movement. Each pair of objects has a 
threshold value specified. If the Euclidean distance between 
this pair of objectss is within the threshold, then they are 
connected by a link. As soon as the pairwise distance goes 
beyond the threshold, the link is broken. The tool was designed 
using OpenGL and C++. OpenGL is a 3D graphics API that 
works with C++ that provides dynamic interaction with the 
user. Mostly it is used for game programming and creating 
3D scenes. For this program, we utilized some of the basic 
features of the API to create an interactive application to 
control the variables of each particular orbiting objects. One 
of the features of OpenGL is the ability to alter the camera 
view. This allows us to see every angle of the orbiting objects 
and rotate around the scene. A snapshot of the visualization 
tool with three moving points is shown in Fig. 19. 

The floor is based on a 24 x 24 grid. The center of each orbit 
has the ability to be moved along the X axis, —12 < X < 12, 
and the Z axis — 12 < Z < 12. The F-axis, which represents 
the height above the floor, allows the object to increase in 
height in the range 1 < F < 10. The radius of each circular 
orbit can be modified in the range 1 < R < 10. Each object 
has a connectivity threshold to each other. Each threshold, 
A-B, B-C,SindC-A has the following range 1 < T < 20. 
Using the standard distance equation between two points in 3D 
space, a link will appear to declare if the objects are within 
the given threshold limit. The speed of each object is based 
on the system clock, which will vary between each computer. 
The speed of each object can be increased up to 5 units. Prior 
to each object being set in motion, they can be positioned 
strategically around their own orbit. The motion of all objects 
can be paused to analyze a specific situation. 

X. Experimental Results and Discussion 

In this section we present the experimental evaluation results 
of two strategies proposed in Section VII. The goal of these 
experiments were to find the impact of change of (i) radius 
of the coverage sphere (rg) and (ii) height of the air corridor 
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(a) Lac = 100, Wac = 70, Hac = 10 (b) Lac = 100, Wac = 180, Hac = 10 (c) Lac = 100, Wac = 100, Hac = 10: Strategy 

1 becomes the same as Strategy 2 



Fig. 20. Variation of the Tq, n and the objective function for Strategy 1 and 2 (all values on y axis) with variable Ts and fixed value of Lac, Wac and Hac 




5 10 15 20 25 30 35 % 10 15 20 25 30 35 5 10 15 20 25 30 35 

(a) Lac = 100, Wac = 70, rs = 20 (b) Lac = 100, Wac = 180, rs = 20 (c) Lac = 100, Wac = 100, Ts = 20: Strategy 1 

becomes the same as Strategy 2 



Fig. 21. Variation of the Tq, n and the objective function for Strategy 1 and 2 (all values on y axis) with variable Hac and fixed value of Lac, Wac and Ts 




Fig. 22. The change of objective function and its components with n and ro 



{Hac) on (a) radius of the circular orbit of the flying ANPs 
(ro), (b) the number of ANPs in each orbit (n), and (c) the total 
number of ANPs (mn) needed to provide complete coverage 
for the air corridor, specified by its length, width and height 
parameters Lac-,Wac-,Hac^ respectively. Fig. 20 show impact 
of changing Vs on Tq, n and mn for different sets of values for 
Laci Waci Hac ^nd for two different strategies 1 and 2. Fig. 21 
show impact of changing Hac on ro, n and mn for different 
set of values for L^c, Wac, and for two different strategies 
1 and 2. The parameter values for L^c, Wac^ Hac and used 
for the experimentation are indicated in the figures. 

Since the optimal coverage problem turned out to be a 
non-linear optimization problem (equations (20) to (23)), we 
used the non-linear constrained program solver Nimbus [19] to 
solve it using two different ANP orbit placement strategies 1 
and 2. In the following we discuss some experimental results, 
some of which are intuitive, some others are not. 

Observation 1: From Fig. 20, it can be seen that increase in 
Vs results in increase in Vq and decrease in mn for both the 
strategies 1 . This is somewhat intuitive as it is only natural to 



expect that as the radius of the coverage sphere increases, the 
radius of the circular orbit of the ANPs will increase and the 
total number of ANPs needed to cover the entire air corridor 
will decrease. It may also be noted that when is too small 
compared to Hac, there may not be a feasible solution. 

Observation 2: From Fig. 21, it can be seen that increase in 
Hac results in decrease in Vq and increase in mn for both 
the strategies 1 and 2. This is also somewhat intuitive, as the 
height of the air corridor increases, the radius of the circular 
orbit of the ANPs has to decrease (please see discussion in 
Section VII) and the total number of ANPs needed to cover 
the entire air corridor must increase. 

Observation 3: From Figs. 20 and 21, it can be seen that, n, the 
number of ANPs in an orbit remains a constant irrespective of 
changes in L^c, Wac^ Hac and rg. This result is not at all obvi- 
ous. However, on closer examination of the objective function 
in equations (20) and (22), one can find an explanation for 
this phenomenon. The plot of n/co5^(^) versus n is shown 
in Fig. 22(a). This factor is present in the objective function 
for both the strategies. From Fig. 22(a), n/co5^(^) reaches its 
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minimum value when n = 5. Therefore the objective functions 
in equations (20) and (22) are nuninuzed when n = 5. This 
nature of n also explains the fact in observation 1, where 
mn decreases when m decreases (i.e., when Tq increases). 
Similarly, it also explains the fact in observation 2, where mn 
increases when m increases (i.e., when decreases). 

Observation 4: From the Figs. 20(a), 20(b), 21(a), 21(b), it can 
be seen that the cost of the solution (i.e., the number of ANPs 
needed to provide complete coverage of the air corridor) using 
strategy 1 is less than that of strategy 2. Although, the reason 
for this phenomenon may not be obvious at a first glance, on 
closer examination, we can explain the phenomenon. Given the 
fact that L^^ + W^^ > 2LacWac and presence of these two 
terms in objective functions of strategies 1 and 2 (equations 
(20) and (22) in page 8), it is not surprising that cost of the 
solution strategy 1 is less than that of strategy 2. 

From our experiments we learn that (i)strategy 1 performs 
better than strategy 2 in all cases, except where Lac = Wac, 
for which both the strategies are identical, (ii)the number of 
ANPs in an orbit remains a constant (5) irrespective of the 
values of Lac^Wac^Hac and r^, when the objective function 
is specified by equations (20) or (22), and (iii)to optimize the 
objective function, the radius of the circular orbit of the ANPs 
should be made as large as possible subject to the constraint 
that the height of the corresponding invariant coverage cylinder 
is at least as large as the height of the air corridor Hac- 

XI. Conclusion 

Existence of sufficient control over the movement pattern 
of the mobile platforms in Airborne Networks opens the 
avenue for designing topologically stable hybrid networks. In 
this paper, we discussed the system model and architecture 
for Airborne Networks (AN). We studied the problem of 
maintaining the connectivity in the underlying dynamic graphs 
of airborne networks with control over the mobility parameters 
and developed an algorithm to solve the problem. 

References 

[1] J. L. Burbank, P. H. Chimento, B. K. Haberman, and W. T. Kasch, "Key 
Challenges of Military Tactical Networking and the Elusive Promise 



of MANET Technology," IEEE Communication Magazine, November 
2006. 

[2] M. Conti and S. Giardano, "Multihop ad-hoc Networking: the Reality," 

IEEE Communications Magazine, April 2007. 
[3] S. M. Alam and Z. J. Haas, "Coverage and Connectivity in Three- 

Dimensional Underwater Sensor Networks," Wireless Communications 

and Mobile Computing, vol. 8, no. 995-1009, 2008. 
[4] , "Coverage and Connectivity in Three-Dimensional Networks," in 

International Conference on Mobile Computing and Networking, 2006. 
[5] H. Tezcan, E. Cayirci, and V. Coskun, "A Distributed Scheme for 3D 

space Coverage in Tactical Underwater Sensor Networks," in MILCOM, 

2004. 

[6] S. Poduri, S. Pattern, B. Krishnamachari, and G. S. Sukhatme, "Sensor 
Network Configuration and the Curse of Dimensionality," in EmNets, 
2006. 

[7] F. Chen, P. Jiang, and A. Xue, "An Algorithm of Coverage Control for 
Wireless Sensor Networks in 3D Underwater Surveillance Systems," 
Advanced Intelligent Computing Theories and Applications, vol. 5226, 
no. 1206-1213, 2008. 

[8] R. Lei, L. Wenyu, and G. Peng, "A Coverage Algorithm for Three- 
Dimensional Large-Scale Sensor Network," in Intelligent Signal Pro- 
cessing and Communication Systems, 2007. 

[9] C. Huang, Y. Tseng, and L. Lo, "The Coverage Problem in Three- 
Dimensional Wireless Sensor Network," in GLOBECOM 2004, 2004. 
[10] "City-wide Wi-fi Projects." [Onhne]. Available: 
http://www.seattlewireless.net,http://www. wirelessphiladelphia.org, 
http://www.waztempe.com 
[11] A. Tiwari, A. Ganguli, and A. Sampath, "Towards a Mission Planning 
Toolbox for Airborne Networks: Optimizing Ground Coverage Under 
Connectivity Constraints," in IEEE Aerospace Conference, March 2008, 
pp. 1-9. 

[12] B. Epstein and V. Mehta, "Free Space Optical Communications Routing 
Performance in Highly Dynamic Airspace Environments," in IEEE 
Aerospace Conference Proceedings, 2004. 

[13] "Doppler Effect." [Onhne]. Available: 

http://en.wikipedia.org/wiki/Doppler_effect 

[14] M. Mesbahi, "On State-dependent Dynamic Graphs and their Control- 
lability Properties," IEEE Transactions on Automatic control, vol. 50, 
no. 3, pp. 2473-2478, 2005. 

[15] , "Controlling Connectivity of Dynamic Graphs," in 44th IEEE 

Conference on Decision and Control, December 2005. 

[16] M. M. Zavlanos and G. J. Papas, "Potential Fields for Maintaining 
Connectivity of Mobile Networks," IEEE Transactions on Robotics, 
vol. 23, no. 4, August 2007. 

[17] T. H. Cormen, C. E. Leiserson, R. L. Rivest, and C. Stein, Introduction 
to Algorithms. McGraw Hill, 2001. 

[18] R. Diestel, Graph Theory. Springer, 2005. 

[19] "Nimbus: Interactive multi-objective optimization system." [Online]. 
Available : http : //w w wnimbus . it. j yu . fi/N4/index . html 



